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Abstract: 

We apply our solution of the nonlinear evolution equation to the case of J/ip photo and DIS 
production on nucleons and nuclei targets. The uncertainty in the J/ij) wave function normal- 
ization due to Fermi motion is treated as a free parameter. We obtain good reproduction of the 
HERA experimental data on a proton target. Calculations of J/ip mesons coherent production 
on nuclei targets are presented and discussed. Our analysis supports the conclusions reached in 
our previous studies, stressing the importance of nonlinear evolution in the kinematical domain 
of high-density QCD. 
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1 Introduction 



It was recognized long ago [1] that the interactions of virtual photons at high-energies and low 
virtualities are governed by nonlinear evolution. Over the years, different theoretical methods 
of calculating processes in the kinematical region of high-density QCD (hdQCD) have been 
proposed 13 IH 13 EI , leading to the construction of an equation, which incorporates both 
linear evolution due to parton splitting, and nonlinear evolution due to the recombination of 
partons at high-density This equation is given in Eq. below, using the notations of [7j, for 
N(r±, x; b), the imaginary part of the amplitude of a dipole of size r±, which scatters elastically 
at impact parameter b. 

In this paper we extend our recent investigations [HI E] of the nonlinear evolution to vector 
mesons, focusing on J/ip photo [TU] and DIS [TT] production. As this process is characterized 
by a scale of the order of the charmed quark mass, m c , it is considered a primary candidate 
for investigating the kinematical region of QCD on the boundary between perturbative and 
hdQCD. 

In our previous publications ^21 El we demonstrated that the cross section experimental 
data for J/ ip production is well reproduced if one employs two damping factors, calculated using 
a Glauber like shadowing-correction (SC) formalism ^3]. Broadly speaking, each damping 
factor may be considered as one iteration in an iterative procedure for determining N(r±, x\ b). 
These iterations correspond to Glauber rescatterings where the first iteration is attributed to 
SC due to the passage of a dipole through the target, and the second iteration is attributed 
to one small x gluon emitted by the dipole prior to the interaction. We believe that the SC 
formalism was, at the time, an important preliminary step towards the more rigorous treatment 
of 0. 

A particular outcome of > is that the J / ip experimental data for the forward differential 
cross section slope, B, can be well reproduced by assuming that the profile function in the 
impact parameter space is the Fourier transform of an electromagnetic dipole-like form factor. 
Our numerical calculations of B produced a satisfactory fit to the experimental data with a 
hadron radius R 2 = 10 GeV~ 2 . On the other hand, in 0, a good fit to the F2 data was obtained 
with R 2 = 3.1 GeV~ 2 for an exponential form factor, and R 2 = 4.5 GeV -2 , for a dipole-like 
form factor. We consider the radii in as effective radii which are reduced due to (i) the 
presence of inelastic diffractive processes; and (ii) the observation that in the simple expression 
typically used for Odipoie, the dipole-proton cross section, a small anomalous dimension has been 
assumed. As we shall see, when calculating the cross section for J/ip production, it is necessary 
to include an effective 6-dependence, which we choose to extract from the HERA experimental 
data. 

This paper is organized as follows: in section El we briefly review the nonlinear equation 
and the approximate numerical solution to it, as obtained in 0; in section El we calculate 
the integrated cross section for J/ip photo- and DIS- production, and compare to the relevant 
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experimental data; in section 0] we present our predictions for the production of J / if) from 
scattering of electrons on heavy nuclei. Our summary and conclusions are given in section |3J 



2 The Dipole Scattering Amplitude 



The imaginary part of the dipole scattering amplitude, N, is a solution of a nonlinear evolution 
equation, which characterizes the low x behavior of the parton densities, while taking into 
account hdQCD effects, thereby obeying the unitarity constraints. The equation describes the 
interaction with a target of a parent dipole, of size x i, and of two dipoles, of sizes x 12 and x 2, 
which were produced by the dipole of size x i. The probability for the decay of the dipole of 
size xqi is given by the square of its wave function, which, in a simplified form, can be written 



as x 



01 



X 



12- 



Each of the produced dipoles can interact with the target independently, with respective 
amplitudes of A r (xi 2 ,?/; b — ^x 2) and A r (x 2,?/; b — ^xi 2 ), where y denotes the rapidity vari- 
able and h the impact parameter. However, adding these contributions clearly overestimates 
the dipole-nucleon interaction, since one must consider the probability that during the inter- 
action, one dipole is in the shadow of the other. This negative correction factor is given by 
-iV(xi2, y;b- ±x 02 )iV(xo2, y; b - ±x 12 ). 

Thus, N can be written in the form 



iV(x i, Y; b) = iV(xoi, Y ; b) exp 



Y 



dyexp 
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2 iV(x 02 , y; b - ^12) - iV"(x 2, y; b - -x 12 )iV(x 1 2, y\ b 



;X 2, 



where, Y = — lnx, Yq = — Itlxq and p is an ultraviolet cutoff, which does not appear in the 
physical quantities. 

The linear part of (^J) is the LO BFKL equation , which describes the evolution of the 
multiplicity of the fixed size color dipoles as a function of the energy Y. The equation sums the 
high twist contributions. Note, that the linear part of also has higher twist contributions 
and the main contribution of the nonlinear part is to the leading twist (see [2] for general 
arguments and JH] for explicit calculations). 

For completeness, we provide a brief description of the steps taken in |H| for obtaining an 
approximate solution for iV(xoi, Y; b). 
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The initial conditions for (0) were taken at xq = 10 2 in the eikonal approximation, ac- 
counting for multiple dipole-target interactions: 

iV(x 01 ,X ;6) = 1 - e -5-input(*Ol,xo)S(b) ; (2) 

where 

2 

a in put(xoi,xo) = ^j^x 2 01 xG BGLAP (x ,A/x 2 01 ), (3) 

and S(b) is the profile function in impact parameter space. As stated, we found that the 
experimental data of the differential cross section slope, B, is well described by the following 
profile function, which is the Fourier transform of electromagnetic dipole-like form factor: 

As a first step, the 6-dependence of ((TJ) was neglected. Thus, all twist contributions for the 
evolution were summed, using the initial condition (J2J) at b = 0. Then, once an approximate 
solution was obtained, the b dependence is restored, assuming similar 6-dependences for both 
the solution and the initial conditions. Specifically, the following ansatz was used for the 
6-dependence of N: 

N( r± ,x;b) = (1 - e -^ ± ) S (b)/s(o)^ ^ 

where 

K(x,r ± ) = - ln(l - N(r±,x,b = 0)). (6) 



3 Cross section for J ftp production 



The cross section for J /if) production is given by: 

a(j*p — ► Vp) = J d 2 b J dz d 2 r±^ 7 *(r ± , z,Q 2 )A(r ± , x;b)^ v (r ± , z) 



(7) 



where A(r±, x; b) is the imaginary part of the production amplitude in impact parameter space, 
and and are, respectively, the wave functions of the J /if) and the virtual photon. 
Formally speaking tyy should depend on b. However, as the dominant contribition of the J /if) 
wave function comes from short distances, this dependence is neglected. The evaluation of ((7j) is 
done by first performing the polarization summation of the J/^-photon overlap function, \Ev x 
\I/ 7 » . This overlap function can be derived in the r± representation using the spin structures of 
the vector meson wave function J7j and the well-known photon wave function |18j . 
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A detailed analysis of the overlap function has already been made in ^J], where it was found 
that, in the momentum represent aion, the overlap function for transverse (T) and longitudinal 
(L) polarized photon is given by [P9*] : 

j / 1 \ j /j \ 2(z 2 + (l-z) 2 )kla 2 + m 2 c (a 2 -k 2 ± ) T . 

V v (k ± , z) x ^ 7 *,t(^±, z) oc -i i ' \ —^v(k ± , z) , (8) 

7Tl c yQ/ ~T~ 

-r /, -r 2z(l — z)Q(a 2 — k 2 , ) T . , n . 

* V (k ± , Z) X %*, L (k ± , Z) OC ^ ( 2 ^2N3 *v(k±, Z) , (9) 

with a 2 = z(l — z)Q 2 + m 2 . Here x denotes the polarization and helicity summation. For a 
given J / ip spatial distribution, the transformation to configuration space is straightforward and 
is given in Eqs. (|10 p -([ll| ) below. Although many models exist for the spatial distributions of 
vector mesons (see for example, [12112111121112211231211), we choose to approximate \&y(rj_, z) to 
be ^v{ r ±. = 0,z = |). As in our previous publications [121 CBj, we consider below a deviation 
from this approximation due to relativistic effects, produced by the Fermi motion of the bound 
quarks within the vector meson. 

The effect of this motion, however, strongly depends on the charmed quark mass, m c . If one 
assumes, for example, that m c = ~ 1.55 GeV, then, by definition, there is no correction 

due to Fermi motion. On the other hand, in [20], it has been assumed that m c ~ 1.50 GeV, 
and a suppression factor of the cross section of about 0.25 was obtained, with almost no energy 
dependence. Hence, as it stands, the contribution of this effect lies within a substantial range 
of uncertainty in which m c varies by no more than 0.05 GeV. 

We, therefore, consider the effect of Fermi motion as an uncertainty of the wave function 
normalization, introducing it as an overall (energy- independent) suppression factor, Kp, which 
we use as a fitting parameter. Using the above approximation we Fourier transform Eqs. (|5|) 
and (El) and obtain: 



£ (c-MO- j*2(C))+m.(^tf 2 (<)-C*,(o)} 



»v(r, = M =i) x * 7 , l(ri ; Q >) = i ^^±§ (^.(O-CJWO). (") 

where ( = ar±, Ki,i = 1, 2 are the modified Bessel functions and T ee = 5.26 KeV is the leptonic 
width of the J/ip. 

Before presenting our results for the cross section ( | 7 | ). we discuss the impact parameter 
dependence of the production amplitude, A. At first sight, A(r±, x; b) equals N(r±, x; b) where 
its 6-dependence is related to the solution of Eq. ((H) at b = through the ansatz of Eq. (JHJ. The 
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Figure 1: Rescatterings for total DIS cross section. 

parameter R in S(b) [see Eq. (0J)] is related to the hadron size and has already been determined 
in [Oj by fitting N to the experimental data of F^. However, as stated in the introduction, a 
good fit to F 2 was obtained using a relatively small value of R 2 = 4.5 GeV 2 , whereas we know 
that the the measured J/ip differential slope, B, is consistent with a radius which is more than 
twice larger. 

Hence, it turns out that A differs from N. The correct 6-dependence of A is given by the 
following expression: 



A(r ± , x;b) = J d 2 b'N(r ± , x; b') S'(b - V) . (12) 

To understand Eq. (|T2|) we need to recall the procedure for calculating the total cross section 
for the deep inelastic process which is shown in Fig.^ If one uses the specific model in which 
the proton consists of two colour dipoles (see Refs. (251 EE] for details), the total deep inelastic 
cross section is given solely by the scattering amplitude of one particular dipole present in 
the proton target. More specifically, as there is no momentum transfer which is involved in 
the process, the total cross section does not depend on the probability of finding a second 
(spectator) dipole in the target. 

On the other hand, in the process of J/i/j production, the momentum transfer, q, is not zero 
(see Fig. EJ). In this process, for fixed momentum transfer, the t = —q 2 dependence of the J/ip 
production relates both to N(r±,x; q) (the amplitude for the photon to scatter off one dipole) 
and to the probability to find the second dipole having a momentum q inside of the recoiled 
proton. Denoting this probability by S'(q), the production amplitude is thus proportional to 
N(r±,x; q)S'(q). A product in momentum representation is equivalent to a convolution in the 
impact parameter, hence the form of Eq. (|12|) . 

Actually, we can estimate S'(q) using the Born approximation for the scatttering amplitude 
shown in the diagram of Fig. El Assuming a simple factorized form of the proton wave function, 

*proton(r±,p±) = 4>(r ± ) 0(p ± ) (13) 
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Figure 2: Rescatterings for J/ip cross section at fixed t = —q 2 



where p± = r\ t ± — r 2i j_ and r±_ = r 3i j_ — (r x j_ + r 2i j_)/2, we can easily calculate S'(q) 

S'(q) = fd 2 r ± \<f>(r ± )\ 2 e^ 1 (14) 



Comparing Eq. ()14j) with the expression for the electromagnetic form factor for proton 



F(q) 



d 2 r ± \(f>(r± i 



12 a 



(15) 



one can see that S'(b) can be described by Eq. (jlj) with R = R proton /2, where -R pro ton is the 
electromagnetic radius of proton. 

Alternatively, we can choose a different strategy and extract the typical size for the profile 
function S'(b) using the experimental data. Generally, the differential slope is related to the 
average of the square of impact parameter, weighted by the the dipole-proton cross section. 
The correct averaging procedure would be to integrate both over the impact parameter and 
over the dipole size, thereby reconstructing the energy dependence of B (the shrinkage of the 
diffraction peak): 

Jd 2 r ± b 2 d 2 bN(r ± ,x,b) 



B x) 



1 IV? 
2 



(16) 



2Jd 2 r ± d 2 bN(r ± ,x,b) 

The value of B as obtained from ()16j) is, of course, not universal for all processes, and it directly 
depends on the choice of S(b). For the purpose of calculating the effective radius to be used in 
S'(b), we define the deviation of \{b 2 ) from the measured slope to be: 



B' = B, 



cxp 



(17) 
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Figure 3: The amplitude in the Born approximation for the J/ip cross section at fixed t = 

„2 



where, -B ex p is taken from experimental data (see, e.g., |2Zj)- S'(b) is then calculated from 
Eq. flU) with the substitution R — ► R' = 2B'. 

The resulting 6-dependence of A(x; b) = J d 2 r±A(r±, x; b)tyj(r±)^v are shown in Fig.0] 
The convolution of S'(b) with N(r±, x; b') is realized by a shift in the maximum of the scattering 
amplitude in the impact parameter space. The decrease of the amplitude in the low b region 
may be understood as a signature of shadowing corrections needed to insure unitarity. 

Our calculations of cr(7*p — > J/ip) were basically in accordance with Eq. (J2J), with modi- 
fications due to the contribution from the real part of the production amplitude [2E] and the 
skewed (off diagonal) gluon distribution [2H|. The contribution of the real part is given by: 

C 2 R =(1+P 2 ), (18) 

where, 

p = ReA/ImA = tg(—) (19) 

2 

and, in our approach, 

A = <91n(iV)/<91n(-). (20) 

x 

Note that with this definition of A, its value obeys the unitarity constraint according to which, 
at a fixed value of Q 2 , A is a decreasing function of the energy with A ^ as i - > j^j. 

The off diagonal contribution is given by: 



ciJ^±m 2 . (21) 
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Figure 4: The impact parameter dependence of the J ftp production amplitude. 
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Formally, the definition (|2U|) of A is r^-dependent. In practice, the corrections due the r±- 
dependence of A are rather small. Thus, A is computed at fixed scale r\ = ^/{Q 2 + M 2 ^). 

The results of our calculations of Eq. (JJJ), given the above modifications, are shown in Fig. El 
for J/ip photoproduction, where we define x = Mj^/W 2 . We found that the HI data [3U| 
are compatible with K F = 0.6 (x 2 /n.d.f. = 0.5) while the ZEUS data require K F = 0.76 
(x 2 /n.d.f. = 0.6). As stated, the factor K F is the deviation of our approximation for the J/ip 
static potential wave function from a more realistic model, in which the spatial distribution of 
J/ip depends on r± and z. Our results are consistent with [21], where a comparison between a 
realsitic and static J/ip wave functions was made. 

The values of K F which were obtained for Q 2 = 0, were used to predict the cross section for 
Q 2 > 0. In Fig.Elwe compare the Q 2 > predictions to the preliminary data which were read 
off plots in [llj. Within the experimental errors, our reproduction of the data is good, except 
for data at Q 2 = 16 GeV 2 where we underestimate the HI point (W ~ 80 GeV) and one of the 
ZEUS points (W ~ 180 GeV). In spite of this, the overall x 2 / n -d.f. for photoproduction and 
DIS data (including these two problematic points) is 0.73 for ZEUS and 0.84 for HI. 



4 Predictions for Heavy Nuclei 

In this section we assess the unitarity effects in J/ip coherent production in photo and DIS on a 
nuclear target. In such processes, unitarity effects are more pronounced than in e — p collisions. 
Thus, with nuclear targets one can access the region of hdQCD at values of x which are larger 
than those characterizing this region at HERA experiments. 

The production of vector mesons on a nuclear target can be calculated in a straightforward 
manner using methods similar to the above. In principle, (0) can be solved separately for each 
atomic number A, with the same form of initial conditions. To this end we choose to present 
our predictions using the Glauber approach, and for the time being to postpone the solution 
of JTJ) for nuclei. 

In the Glauber approach, the cross-section for the J / ip production, including contributions 
for the percolation through a nucleus of both a \qq) system and a \qqg) system, is given by the 
following expression: 

J dz d 2 r ± ^ r (r ± , z,Q 2 )^ J / i ,{r±, z) x 



a(j*A — ► J/ijA) = / d 2 b 



(22) 
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Figure 5: The ZEUS and HI J/ip photoproduction data compared with our calculations 
using NLE. 
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Figure 6: J/ip DIS and photoproduction, data and our calculations using NLE. The upper 
(lower) curves at each Q 2 bin correspond to different Kp values, which were fitted 
to ZEUS (HI) photoproduction data. 
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where the opacities Vt q and Q g are defined as: 



Q q (b,r±,x) = ^r\a B xG (x, S A (b) , (23) 

9 

Q g (b,r ± ,x) = -Sl q (b,r±,x), (24) 

Sa(&) is the number of nucleons in a nucleus interacting with the incoming dipole and xG is 
the gluon density obtained from linear DGLAP evolution (we used CTEQ parameterization 
pTT] for this particular computation). 

We are cognisant of the fact that when dealing with heavy nuclei interactions, neither a 
dipole nor a Gaussian are adequate for 5^(6). When the probing hadronic system percolates 
through a heavy nucleus, it experiences multiple rescatterings over a substantially long region 
of the impact parameter space. On the other hand, once the system leaves the heavy nuclei, 
the number of strong interactions decreases rapidly. 



Our approach is to use the Wood-Saxon parameterization [32J for the impact parameter 
dependence, 5^(6): 

S A (b) = P [ dr l RA (25) 



1 + e—tr 

where rn is the longitudinal distance from the target, h and Ra are parameters which are taken 



from experimental tables [323 > r = y r f\ + b 2 and the normalization factor p is defined in the 
following relation: 

r dr\\ d 2 b . , 

P / W = ^- (26) 

J 1 + e^ 

Note that Ra is associated with the radius of a nucleus with atomic number A. 

To determine the variance between nuclei targets and a nucleon target, we need to examine 
the A dependence of the integrated cross section, by comparing the calculated cross section to 
a power behaviour of a oc A s . The maximal value of the exponent, S max , depends on the profile 
function 5*^(6), and can be found from the following integral: 



A 5 " 



Jd%\S A {b)\ 2 (27) 



For an exponential profile, 5 max = | and for a Wood-Saxon profile of (|23|). 5 max ~ 1.43. A 
deviation from the maximal value of 5 would be a clear signature of saturation of the growth 
of the cross section. Fig.[7^i shows the cross section as a function of A with an arbitrary 
normalization, for x = 10 -5 , 10 -4 and 10~ 3 . To illustrate the saturation effect, we compare the 
curves to a power behavior of a oc A^ 3 . Unitarity conserving effects are already appreciable at 
x = 1(T 3 . 
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Eq. (|22|) takes into account the interaction of the quark- ant iquark pair and the fastest gluon. 
On the other hand, we can improve our Glauber approach by using the solution of the non- 
linear equation, since the solution to the non-linear equation takes into account all possible 
interactions of partons with the nucleon target. As stated, R 2 , the effective radius of nucleon 
turns out, in our approach, to be very small. With such a small radius, the typical parameter 
that governs the dipole interaction with the nuclear target is also very small, namely, 

R 2 p S A (0) < 1. (28) 

In the standard approach, however, R 2 5^(0) was considered to be large. We thus have to take 
into account all rescatterings inside the proton in the first stage of our approach and treat the 
interactions with different nucleons in a nucleus using the Glauber formula. 

Specifically, the cross section can be written in the form: 
a(^*A — ► J/ipA) = 



J d 2 b J dzd 2 r ± V r (r ± ,z,Q 2 )y J/i ,(r ± ,z) (l - e 



i<7 (dipole— nucleon) Sa{V) 



(29) 



where a(dipole — nucleon) is equal to 2 / d ?bN(r , x\ b) where A" is dipole-nucleon amplitude. 

The following figures illustrate the effect of the inclusion of all rescatterings. Fig.[7)D shows 
the cross section of Eq. ()29|) as a function of A with an arbitrary normalization, for x = 
10" 5 , 10" 4 and KT 3 . 

In Fig. [7k we display the characteristic exponent, 5, as a function of x, for the calculated 
cross-sections of Eqs. (l22|) and (J29|) . The upper curve in Fig. corresponds to the numerical 
calculation of Eq. (|29jl and the lower curve corresponds to Eq. ()22jl . In terms of S, the differences 
between the two approaches are negligible. On the other hand, as further detailed below, there 
is a considerable normalization difference. 

To investigate the normalization difference, we define the following ratio: 



* = (3°) 

Eq. *™ 



We have calculated 1Z for J/ip photoproduction on a Gold target (A ~ 200). The results, as 
a function of the energy, for several values of Q 2 , are shown in Fig.UJi. Our calculations show 
that the full calculation is suppressed by a factor of about 0.4 with respect to the calculation 
in which only a \qq) and a \qqg) systems are considered to interact with the nucleus. 

Our prediction for the energy dependence for J/ip production on a Gold target, in the 
calculation which includes all parton interactions is shown in Fig. |H] for several values of Q 2 . In 
this figure, we have used the wave function normalization which we have extracted from the 
data on proton target. Due to the uncertainty of this normalization, and to the 0.4 ratio of 
Eq. (jnnjL our rough estimation of the relative errors of the curves are of the order of 15-20% . 
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Figure 7: J/ip photoproduction on a nuclear target: (a)-(b) show the ^4-dependence of the 
cross section for different values of x, normalized to unity at A = 50 [(a) shows 
the contribution due to the interactions of \qq) and \qq~g) with the nucleus, and 
(b) shows the contributions of all interactions]; (c) shows the values of 5, as a 
function of x, where the lower curve corresponds to plot (a) and the upper curve 
corresponds to plot (b); and (d) shows the normalization ratio between plot (a) 
and plot (b) as explained in the text. 
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Figure 8: Predictions fof J /if) photo- and DIS production on a Gold target, using the 
Glauber approach of Eq. (|29jl for different values of Q 2 as a function of W. 
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5 Summary 



We have extended our recent investigation of the approximate solution of the nonlinear evolution 
equation to the case of J / ip photo- and DIS production. We have convoluted our previous ansatz 
for the 6-dependence of the amplitude with an additional profile which is extracted from the 
electromagnetic form factor. The resulting impact parameter dependence of the amplitude 
exhibits a decrease near 6 = 0. We believe that this decrease is due to the deviation from the 
linear evolution equations, and is a signature for the onset of unitarity taming effects. 

We have used the solution to the NLE and obtained a satisfactory reproduction of both the 
photo- and the DIS production data. Our free parameter was the uncertainty in the J j ip wave 
function normalization due to Fermi motion of the cc pair. In addition, we took into account 
modifications of the cross section due to the real part of the amplitude and the existence of 
skewed gluons within the nucleon. 

We used the Glauber approach to calculate the cross section for J / ip production on a nuclear 
target. Our predictions demonstrate that such processes are important for the understanding 
of unitarity taming effects. We have found that these effects start to dominate for x < 10~ 3 . 
Hence, the predictions presented here may be verified with the future eRHIC experiments. 

It is important to note that, apart from the uncertainities due to vector meson wave func- 
tion, the calculations on nuclei do not contain any fitting parameters. More specifically the 
uncertainity due to b dependence of the solution for nuclei is quite small. In other words, the 
initial conditions for the BK equation jHJ U\ are determined solely by the proton DIS data. 
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